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We discuss the production of top-anti-top quark pairs in association with a hard jet at the Tevatron 
and at the LHC and we report on the calculation of the next-to-leading order QCD corrections to 
this process. Numerical results for the tt-njet cross section and the forward-backward charge 
asymmetry are presented. The corrections stabilize the leading-order prediction for the cross 
section. In contrast, the charge asymmetry receives large corrections. The dependence of the 
cross section as well as the asymmetry on the minimum transverse momenta used to define the 
additional jet is studied in detail for the Tevatron. 
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1. Introduction 

The top quark is the heaviest of the known elementary particles. More than ten years after 
its discovery, the dynamics and many properties of the top quark, such as its electroweak quantum 
numbers, are not yet precisely measured. It is widely believed that the top quark plays a key role 
in extensions of the Standard Model. This renders experimental investigations of the top quark 
particularly important. Up to now the main (direct) source of information on top quarks are top- 
quark pairs produced at the Tevatron. Only recently first evidence for single-top production has 
been found [|I]]. It is important to note that in the inclusive tt sample a significant fraction comprises 
tt-f-jet events. An investigation of the process of tt production in association with a hard jet can thus 
improve our knowledge about the top quark. 

In this context, the forward-backward charge asymmetry of the top (or anti-top) quark [Q, 
|3|, ^, ||] is of particular interest. In inclusive tt production it appears first at one loop, because it 
results from interferences of C-odd with C-even parts of double-gluon exchange between initial 
and final states. This means that the available prediction for tt production — although of one-loop 
order — describes this asymmetry only at leading-order (LO) accuracy. In tt+jet production the 
asymmetry appears already in LO. Thus, the next-to-leading order (NLO) calculation described in 
the following provides a true NLO prediction for the asymmetry. Our calculation will, therefore, 
be an important tool in the experimental analysis of this observable at the Tevatron where the 
asymmetry is measureable as discussed in Ref. [^. 

Measuring the cross section of the related process of tt+y production provides direct access 
to the electric charge of the top quark. Obviously NLO QCD predictions to this process aie im- 
portant for a reliable analysis. They can be obtained from tt+jet production presented here via 
simple substitutions. Finally, a signature of tt in association with a hard jet represents an impor- 
tant background process for searches at the LHC, such as the search for the Higgs boson in the 
weak-vector-boson fusion or ttH channels. 

The above-mentioned issues clearly underline the case for an NLO calculation for tt+jet pro- 
duction at hadron colliders. We report here on a first calculation of this kind as presented in Ref. [^]. 

2. Details of the NLO calculation 

At LO, hadronic tt+jet production receives contributions from the partonic processes qq — > tig, 
qg -^ tiq, qg — > tt^, and gg -^ tig. The first three channels are related by crossing symmetry to the 
amplitude -^ tiqqg. Evaluating 2 — > 3 particle processes at the NLO level, is non-trivial, both 
in the analytical and numerical parts of the calculation. In order to prove the correctness of our 
results we have evaluated each ingredient twice using independent calculations based — as far as 
possible — on different methods, yielding results in mutual agreement. 

2.1 Virtual corrections 

The virtual corrections modify the partonic processes that are already present at LO. At NLO 
these corrections are induced by self-energy, vertex, box (4-point), and pentagon (5-point) correc- 
tions. The most complicated diagrams are the pentagon diagrams. 
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Version 1 of the virtual corrections is essentially obtained following the method described in 
Ref. [0], where ttH production at hadron colliders was considered. Feynman diagrams and ampli- 
tudes have been generated with the FeynArts package [^, ||] and further processed with in-house 
Mathematica routines, which automatically create an output in Fortran. The IR (soft and coUinear) 
singularities are analytically separated from the finite remainder as described in Refs. |Q, [I^. The 
tensor integrals appearing in the pentagon diagrams are directly reduced to box integrals following 



Ref. []11[]. This method does not introduce inverse Gram determinants in this step, thereby avoid- 
ing notorious numerical instabilities in regions where these determinants become small. Box and 



lower-point integrals are reduced a la Passarino-Veltman [12] to scalar integrals, which are either 
calculated analytically or using the results of Refs. [|l^, 14, 15]. Sufficient numerical stability is 
already achieved in this way. Nevertheless the integral evaluation is currently further refined by 
employing the more sophisticated methods described in Ref. [ [l^ ] in order to numerically stabilize 
the tensor integrals in exceptional phase-space regions. 
Version 2 of the evaluation of loop diagrams starts with the generation of diagrams and ampli- 



tudes via QGRAF []17[], which are then further manipulated with Form [ ]18[ ] and eventually automat- 
ically translated into C++ code. The reduction of the the 5-point tensor integrals to scalar integrals 



is performed with an extension of the method described in Ref. [|19p. In this procedure also in- 
verse Gram determinents of four four-momenta are avoided. The lower-point tensor integrals are 
reduced using an independent implementation of the Passarino-Veltman procedure. The IR-finite 



scalar integrals are evaluated using the FF package [20, 21]. 



2.2 Real corrections 



ttqqq'q' 



The matrix elements for the real corrections are given by ^ iigggg, -^ iiqqgg, 
and -^ tiqqqq. The various partonic processes are obtained from these matrix elements by all 
possible crossings of Ught particles into the initial state. 

The evaluation of the real-emission amplitudes is performed in two independent ways. Both 



evaluations employ the dipole subtraction formalism [g2|, g3|, |24[] for the extraction of IR singular- 
ities and for their combination with the virtual corrections. 
Version 1 results from a fully automated calculation based on helicity amplitudes, as described in 



Ref. [25]. Individual helicity amplitudes are computed with the help of Berends-Giele recurrence 



relations [26]. The evaluation of color factors and the generation of subtraction terms is automated. 
For the channel gg -^ iigg a dedicated soft-insertion routine [ 27 ] is used for the generation of the 
phase space. 

Version 2 uses for the LO 2 — > 3 processes and the gg -^ ttgg process optimized code obtained 
from a Feynman diagramatic approach. As in version 1 standard techniques like color decompo- 
sition and the use of helicity amplitudes are employed. For the 2^4 processes including light 
quarks, Madgraph [^] has been used. The subtraction terms according to Ref. [ 24 ] are obtained in 
a semi-automatized manner based on an in-house library written in C++. 



3. Numerical results 



In the following we consistently use the CTEQ6 ]|2^, 3C ] set of parton distribution functions 
(PDFs). In detail, we take CTEQ6L1 PDFs with a 1-loop running a^ in LO and CTEQ6M PDFs 
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Figure 1: Scale dependence of the LO and NLO cross sections for tt+jet production at the Tevatron (left) 
and at the LHC (right) as taken from Ref. |^], where the renormalization scale (/ir) and the factorization 
scale {}Xf) are set equal to /i. 
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Figure 2: Scale dependence of the LO and NLO forward-backward charge asymmetry of the top quark in 
pp -^ tt+jet+X at the Tevatron as taken from Ref. [g| with jx = Hf = jir- 



with a 2-loop running a^, in NLO. The number of active flavours is Np = 5, and the respective QCD 
parameters are Aj^ = 165 MeV and A^^ = 226 MeV. Note that the top-quark loop in the gluon 
self-energy is subtracted at zero momentum. In this scheme the running of a^ is generated solely 
by the contributions of the light quark and gluon loops. The top-quark mass is renormalized in the 
on-shell scheme, as numerical value we take mt = 174 GeV. 

1 for the definition of the tagged hard jet. Un- 
20 GeV for the hardest 



We apply the jet algorithm of Ref. [31] with R 
less stated otherwise we require a transverse momentum of pj.jet > Pj 
jet. The outgoing (anti-)top quarks are neither affected by the jet algorithm nor by the phase-space 
cut. Note that the LO prediction and the virtual corrections are not influenced by the jet algorithm, 
but the real corrections are. 

In Figure [T] the scale dependence of the NLO cross sections is shown. For comparison, the 
LO results are included as well. As expected, the NLO corrections significantly reduce the scale 
dependence compared to LO. We observe that arround /i « mt the NLO corrections are of moderate 
size for the chosen setup. 

We have also studied the forward-backward charge asymmetry of the top quark at the Teva- 
tron. In LO the asymmetry is defined by 
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aLo(3't>0)±aLo(3't<0), 



(3.1) 



where jt denotes the rapidity of the top quark. Cross-section contributions a(yt < 0) correspond to 
top quarks in the forward or backward hemispheres, respectively, where incoming protons fly into 
the forward direction by definition. Denoting the corresponding NLO contributions to the cross 
sections by Sa^^o' ^^ define the asymmetry at NLO by 
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FB.NLO 



^LO 
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NLO 
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(3.2) 



i.e. via a consistent expansion in a^. Note, however, that the LO cross sections in Eq. ( [3.2| ) are 
evaluated in the NLO setup (PDFs, oQ. 

Figure ^ shows the scale dependence of the asymmetry at LO and NLO. The LO result for the 
asymmetry is of order a^ and does therefore not depend on the renormalization scale. The plot 
for the LO result shows a mild residual dependence on the factorization scale, but the size of this 
variation does not reflect the theoretical uncertainty, which is much larger. The NLO corrections 
to the asymmetry are of order a] and depend on the renormalization scale. It is therefore natural 
to expect a stronger scale dependence of the asymmetry at NLO than at LO, as seen in the plot. 
The size of the asymmetry, which is about —7% at LO, is drastically reduced by the NLO correc- 
tions. To investigate the origin of the large NLO corrections to the asymmetry we have studied the 
dependence on the cut value p'^^ used to define the minimal px of the additional jet. The results 
are shown in Table |l]. We observe that both the NLO cross section as well as the NLO asymmetry 
dependent strongly on p'^^. This is related to the fact that the cross section becomes ill-defined in 
the limit p'^^ —>■ due to the appearance of IR divergencies. On the other hand, the LO prediction 
for the asymmetry shows only a mild dependence on p^^. 





cross section [pb] 


charge asymmetry [%] 
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Table 1: Cross section and forward-backward charge asymmetry at the Tevatron for different values of p^^ 
used to define the minimal transverse momentum pj of the additional jet {jX = iXf = jXr = nit). The upper 
and lower indices are the shifts towards pL = mf/2 and pL = Init. 



4. Conclusions 



Predictions for tt+jet production at hadron colliders have been reviewed at NLO QCD. For the 
cross section the NLO corrections drastically reduce the scale dependence of the LO predictions, 
which is of the order of 100%. The charge asymmtry of the top quarks, which is going to be mea- 
sured at the Tevatron, is significantly decreased at NLO and is almost washed out by the residual 
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scale dependence. In addition we have also studied the p™' -dependence of the NLO predictions. 
Further refinements of the precise definition of the forward-backward asymmetry are required to 
stabilize the asymmetry with respect to higher order corrections. 
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